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^ ; Abstract 
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We study medium modifications of J/ip production in cold nuclear media 
in deuterium-nucleus collisions. We discuss several parameterizations of the 
modifications of the parton densities in the nucleus, known as shadowing, 
[ — . an initial-state effect. We also include absorption of the produced J/ip by 

nucleons, a final-state effect. Both spatially homogeneous and inhomoge- 
neous shadowing and absorption are considered. We use the number of binary 
nucleon-nucleon collisions as a centrality measure. Results are presented for 
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d+Au collisions at ^/Snn = 200 GeV and for d+Pb collisions at y/ Snn = 6.2 
TeV. To contrast the centrality dependence in pA and dA collisions, we also 



present pPb results at y/SjyN = 8-8 TeV. 
Ph. 

Q . 

The nuclear, or A, dependence of J/ip production is an important topic of study for 
nuclear collisions. It is essential that the A dependence be understood in cold nuclear 
matter to set a proper baseline for quarkonium suppression in AA collisions. At fixed target 
energies, NA50 [1] has studied the J/ip A dependence and attributed its behavior to J/ip 
break up by nucleons in the final state, referred to as nuclear absorption. However, it is 
also known that the parton distributions are modified in the nucleus relative to free protons. 
This modification, referred to here as shadowing, is increasingly important at higher energies. 
This paper studies the interplay of shadowing and absorption in deuteron-gold collisions at 
the Relativistic Heavy Ion Collider (RHIC) and in deuteron-lead and proton-lead collisions at 
the future Large Hadron Collider (LHC). We consider both spatially homogeneous (minimum 
bias collisions) and inhomogeneous (fixed impact parameter) results. When possible, we 
discuss the results in the context of data, in particular, the PHENIX data from RHIC. 

The nuclear quark and antiquark distributions have been probed through deep inelastic 
scattering (DIS) of leptons and neutrinos from nuclei. These experiments showed that 
parton densities in free protons are modified when bound in the nucleus [2]. At low parton 
momentum fractions, x < 0.02, the ratio of the nuclear parton density relative to the nucleon 
is less than unity (the shadowing region). In the intermediate x regime, 0.02 < x < 0.1, 
this ratio is larger than unity (antishadowing) while at higher x it drops below unity once 
more (the EMC region). The effect also depends on the scale of the interaction, the square 
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of the momentum transfer, Q 2 . As the collisions energy increases, the values of x probed 
in the collision are decreased, making shadowing more important. Since shadowing affects 
the parton distribution functions before the collision that produces the J /if) occurs, it is an 
initial state effect. 

Shadowing should depend on the spatial position of the interacting parton within the 
nucleus [3]. Although DIS experiments are typically insensitive to this position dependence, 
some spatial inhomogeneity has been observed in vN scattering [4]. Thus the effect should 
be sensitive to the impact parameter, b, at which the collision occurs so that the results 
depend on the collision centrality. Central collisions with low impact parameter should 
exhibit stronger shadowing effects than collisions in the nuclear periphery. In Refs. [3,5], 
we discussed the effects of spatially inhomogeneous shadowing on J /if) production in AA 
and dA collisions respectively. In this paper, we augment the dA predictions by including 
final-state absorption of the J /if) by nucleons. 

At lower fixed-target energies, such as those at the CERN SPS, the x values are rather 
high, x ~ 0.18 at midrapidity, and shadowing effects are small. In this energy and rapidity 
regime, the A dependence of J /if) production is attributed to nuclear absorption of the J /if): 
either the pre-resonant cc pair or the J /if) itself can interact with nucleons along its path 
and break up into charm hadrons in reactions such as pJ /if) — > DDX or pj/ip — > DA C . 
However, at energies such as those available at RHIC and the LHC, especially away from 
midrapidity, much smaller x values may be reached, making shadowing more important 
relative to absorption. 

The cc pair that produces the J /if) can be created in either color singlet or color octet 
states. In Ref. [6], absorption was described in terms of the singlet and octet components 
of the J /if} wavefunction, 

I J If) = aol(cc)i) + ai | (cc) 8 £) + a 2 \{cc) 1 gg) + a' 2 \(cc) s gg) + ■■■ . (1) 

In the color singlet model [7], only the first component is nonzero for direct J /if) production. 
The cc pairs pass through nuclear matter in small color singlet states and reach their bound 
state size outside the nucleus, except at sufficiently negative rapidity where the J /if) may 
still hadronize inside the nuclear medium. If cc pairs are predominantly produced in color 
octet states, then the \ (cc)$g) state interacts with nucleons. The produced color octet cc can 
neutralize its color by a nonperturbative interaction with a gluon. Since the \ (cc)sg) state is 
fragile, a gluon exchange between it and a nucleon would separate the (cc)s from the gluon 
and break up the unbound octet [6]. If the \(cc)sg) state evolves without interaction, such 
as in pp collisions, the additional gluon is absorbed by the octet cc pair, hence 'evaporating' 
the color. The color evaporation model (CEM) thence does not care about the relative 
coefficients in Eq. (1). On the other hand, the NRQCD approach, applied to the total cross 
section in Ref. [8], provides the leading coefficients in the expansion of the wavefunction 
in Eq. (1), encompassing both singlet and octet production and absorption. We calculate 
absorption of color singlets and color octets in the CEM as well as the combination of the 
two in NRQCD. The NRQCD approach fixes the fraction of charmonium states produced 
in color singlets and color octets, determining the rate of singlet and octet absorption in 
Eq. (1). 

The nuclear dependence of hard process production such as heavy quarks and quarkonium 
in AB collisions is typically parameterized as 
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&ab = o-NN(AB) a 
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where A and B can be either protons or nuclei and u^n is the production cross section in 
a nucleon-nucleon collision. At central values of rapidity, y, or Feynman x, xp, where xp = 



^Pz/V^nn — 2mr sinh y / v^Sjvtv , m T — \Jp T + m 2 where pp is the transverse momentum 

and m is the particle mass, the A dependence of open charm production is, to a good 
approximation, linear with a = 1 [9]. However, for quarkonium production, a < 1 at 
midrapidity. The most recent measurements by experiment E866 at Fermilab [10] and 
NA50 at CERN [1] suggest a ~ 0.94 — 0.95. Dependence on the kinematic variables xp and 
Pt has been observed. Typically cx(xf) decreases with increasing xf at xf > 0. 

It is as yet unknown if a is a strong function of energy, as predicted by some absorption 
models [11,12]. If color singlet absorption is at work, the absorption contribution should 
decrease with energy because the singlet state will stay small until far outside the target 
[11]. On the other hand, if the absorption cross section depends on the NJ/ip center of 
mass energy, at higher values of ^/Snn the average NJ/ip center of mass energy increases, 
thus increasing the absorption cross section [12]. Since both initial and final-state effects 
such as shadowing and absorption may be dependent on ^/Snn, empirically it would seem 
that a should be energy dependent. At low energies, most analyses have assumed that 
absorption by nucleons is the only contribution to the A dependence. While this may be a 
good approximation at midrapidity, measured by NA50, away from the central region a is 
xp dependent, as previously mentioned. Thus absorption alone is not enough to explain the 
xp dependence, as already noted a number of times, see e.g. Refs. [13-15]. 

Indeed, the characteristic shape of a{xp) at high xp, xp > 0.25, also cannot be explained 
by shadowing alone [13]. However, at heavy ion colliders, the relationship between xf, 
rapidity, and y/SjyN means that this interesting xf region is pushed to very far forward 
rapidities. At y/S^N — 200 GeV and y = 2.5, the forward edges of the PHENIX muon 
arms at RHIC, xf ~ 0.19 for a J/ip with pp = while at \/Snn = 5.5 TeV and y — 4, 
the forward edge of the ALICE muon arm at the LHC, x F ~ 0.03. Therefore, at collider 
energies, a combination of absorption and shadowing effects may be sufficient to address the 
J /if) data. 

In this paper, we discuss the combined effects of shadowing and absorption both in 
minimum bias dA collisions and as a function of centrality at RHIC and the LHC. We focus 
on d+Au collisions at y/S NN = 200 GeV at RHIC and d+Pb collisions at y/S NN = 6-2 TeV 
at the LHC. While it is unclear whether pA or dA collisions will be used at the LHC, one 
advantage of dA is that the energy is closer to that of Pb+Pb collisions at \/Snn = 5.5 TeV 
whereas the p+Pb center of mass energy per nucleon would be 8.8 TeV. Thus the d+Pb 
combination has been suggested as a baseline measurement at the LHC [16]. 

Our J/ip calculations generally employ the color evaporation model which treats all 
charmonium production identically to cc production below the DD threshold, neglecting 
color and spin [17]. The leading order (LO) rapidity distributions of J/^'s produced in dA 
collisions at impact parameter b is 





(3) 
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The partonic cross sections for LO gluon fusion and qq annihilation are given in Ref. [18], 
M 2 = xix 2 Snn and x lf2 = (M/^/Snn) exp(±y). The fraction of cc pairs below the DD 
threshold that become J/i/j's, Fj/^, is fixed at next-to-leading order (NLO) [17]. Both this 
fraction and the theoretical K factor, K th , which scales the LO cross section to the NLO 
value, drop out of the ratios. The K factor is independent of rapidity [19]. While our 
calculations are LO, the shadowing ratios are independent of the order of the calculation 
[5]. We use m c = 1.2 GeV and Q = M [17] with the MRST LO parton densities. We will 
compare our minimum bias CEM results with calculations employing the NRQCD approach 
and discuss any differences in the resulting shadowing and absorption patterns. 

We assume that the nuclear parton densities, F^(x,Q 2 ,r,z), are the product of the 
nucleon density in the nucleus, pa(s), the nucleon parton density, f^x, Q 2 ), and a shadowing 
ratio, Sp S (A, x, Q 2 , f, z), where r and z are the transverse and longitudinal location of 
the parton in position space. The first subscript, P, refers to the choice of shadowing 
parameterization, while the second, S, refers to the spatial dependence. Most available 
shadowing parameterizations, including the ones used here, ignore effects in deuterium. 
However, we take the proton and neutron numbers of both nuclei into account. Thus, 

F?(x,Q 2 ,r,z)=p d (s)f l N (x,Q 2 ) (4) 
Ff(x, Q 2 , b - f, z') = p A (s')Si tS (A, x, Q 2 , b - f, z')ff (x, Q 2 ) , (5) 

where s = \/r 2 + z 2 and s' = \J \ b — f| 2 + z' 2 . With no nuclear modifications, 
Sp S (A, x, Q 2 , f, z) = I. The nucleon densities of the heavy nucleus are assumed to be 
Woods-Saxon distributions with i?Au = 6.38 fm and i?pb = 6.62 fm [20]. We use the 
Hulthen wave function [21] to calculate the deuteron density distribution. The densities are 
normalized so that / d 2 rdzpA(s) = A. We employ the MRST LO parton densities [22] for 
the free nucleon. 

We have chosen shadowing parameterizations developed by two groups which cover ex- 
tremes of gluon shadowing at low x. The Eskola et al. parameterization, EKS98, is based on 
the GRV LO [23] parton densities. Valence quark shadowing is identical for u and d quarks. 
Likewise, the shadowing of u, d and s quarks are identical at Q\. Shadowing of the heavier 
flavor sea quarks, s and higher, is, however, calculated and evolved separately at Q 2 > Q^. 
The shadowing ratios for each parton type are evolved to LO for 1.5 < Q < 100 GeV and 
are valid for x > 10~ 6 [24,25]. Interpolation in nuclear mass number allows results to be 
obtained for any input A. The parameterizations by Frankfurt, Guzey and Strikman (FGSo, 
the original parameterization, along with FGSh and FGS1 for high and low gluon shadow- 
ing) combine Gribov theory with hard diffraction [26]. They are based on the CTEQ5M 
[27] parton densities and evolve each parton species separately to NLO for 2 < Q < 100 
GeV. Although the x range is 10~ 5 < x < 0.95, the sea quark and gluon ratios are unity 
for x > 0.2. The EKS98 valence quark shadowing ratios are used as input since Gribov 
theory does not predict valence shadowing. The FGSo parameterization is available for four 
different values of A: 16, 40, 110 and 206 while FGSh and FGS1 also include A = 197. We 
use A = 206 for Au with FGSo and A = 197 for FGSh and FGS1. 

Figure 1 compares the four homogeneous ratios, S^ KS and S FGSi for Q = 2m c . The FGS 
calculation predicts more shadowing at small x. The difference is especially large for gluons. 
At very low x, the gluon ratios for FGSo and FGSh are quite similar but in the intermediate 
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x regime, the FGSh parameterization drops off more smoothly. On the other hand, the FGS1 
parameterization levels off at a higher value of S P than the other two FGS parameterizations. 
In the antishadowing regime, FGSh and FGS1 are rather similar to EKS98 while FGSo has a 
larger antishadowing effect. Obviously shadowing alone will give an effective A dependence 
as a function of rapidity with y > corresponding to low x, effectively mirroring the curves 
in Fig. 1. 

Figure 2 shows the average value of X2, the fraction of the nucleon momentum carried by 
the interacting parton that produces the J/ip, for pA collisions at the CERN SPS, y/Sj^ = 
19.4 GeV, and dA collisions at RHIC, = 200 GeV, and the LHC, y/S^ = 6.2 TeV. 

At midrapidity at the SPS, (x 2 ) ~ 0.1, in the antishadowing region. Note that if some 
level of antishadowing is indeed present at this energy, the effective absorption cross section 
determined from minimum bias pA data would actually be underestimated, as discussed in 
the context of Pb+Pb collisions in Ref. [28]. At midrapidity at RHIC, (x 2 ) ~ 0.01, near the 
point where Sp < 1. In the forward region, (x 2 ) ~ 10~ 3 at y ~ 2, clearly in the effective 
low x regime due to gg dominance while for negative rapidity, y ~ —2, (x 2 ) ~ 0.1, in the 
antishadowing region. Even lower values of (x 2 ) are reached in the measurable rapidity 
region of the LHC, at y ~ 4, (x 2 ) ~ 10 -5 , deep in the shadowing region. Systematic studies 
over the widest possible (x 2 ) range can map out the gluon distribution, especially for the 
nucleus, if absorption effects are understood. 

We now turn to the spatial dependence of the shadowing, discussed in detail in Ref. [5] . 
Previously, when no inhomogeneous shadowing parameterizations were available, we con- 
sidered two forms of the spatial dependence, one proportional to the local nuclear density, 
Sp ws , and the other proportional to the parton path through the nucleus, Sp . Here we 
show only results for Sp p and compare these to the results with the inhomogeneous FGS 
parameterizations, FGSh and FGS1. The Sp parameterization is, in general, [5] 

S P<p (A,x,Q 2 ,r,z) = l + N p (S P (A,x,Q 2 ) - l)L^^A ( 6 ) 

where N p is chosen so that (1/A) J d 2 rdzp A (s)S Pp (A, x, Q 2 , r, z) = Sp(A,x,Q 2 ). When 
s = \Jr 2 + z 2 ^> Ra-, the nucleons behave as free particles while the modifications are larger 
than the average value S P (A, x, Q 2 ) in the center of the nucleus,. The normalization requires 
(l/A)Jd 2 rdzp A (s)S P:P = S P . 

While there are three homogeneous FGS parameterizations, no spatial dependence is 
provided for FGSo, the case with the strongest gluon shadowing. Therefore, we also use 
Sp with this parameterization as well as with EKS98. The inhomogeneous FGSh and 
FGS1 parameterizations are defined over 0.587 < s < 10 fm for A = 197 and 208. Below 
s min = 0.587 fm, the ratios are fixed to those at s min while S FGSh { (A, x, Q 2 , s) — > 1 for s > 10 
fm. They do not consider the longitudinal spatial dimension. The normalization requirement 
is similar to that of Eq. (6), (1/A) J d 2 sT A (s)S FGSh l (A,x,Q 2 ,s) = S FGSh l (A,x,Q 2 ). Due to 
this fixed maximum, excluding the tail of the nuclear and deuteron density distributions, 
these parameterizations will have a somewhat stronger spatial dependence in a full dA 
calculation. 

Figure 3 compares S FGSo WS and S FGSop with S FGSh (b) at similar values of the homo- 
geneous shadowing ratios. We see that S FGSo p is quite compatible with the available FGS 
inhomogeneous parameterizations. 
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To implement nuclear absorption on J/ip production in dA collisions, the production 
cross section in Eq. (3) is weighted by the survival probability, S ahs , so that 

S ^(b- s,z') = expj- jT 'dz"p A (b- s,z")a ahs (z"-z')} (7) 

where z' is the longitudinal production point, as in Eq. (5), and z" is the point at which the 
state is absorbed. The nucleon absorption cross section, <r a bs, typically depends on where the 
state is produced in the medium and how far it travels through nuclear matter. If absorption 
alone is active, i.e. Sp s = 1, then an effective minimum bias A dependence is obtained after 
integrating Eqs. (3) and (7) over the spatial coordinates. If S ahs = 1 also, a dA = 2Aa pN . 
When S ahs ^ 1, a dA = 2A a a pN where, if a abs is a constant, independent of the production 
mechanism for a nucleus of p A = Pq6{R A — b), 

a=l-^L (8) 

where r = 1.2 fm. Although we assume absorption is only effective for the heavy nucleus, 
the spatial dependence of the deuteron wavefunction is included in Eq. (3). The impact 
parameter dependence of absorption alone is shown in Fig. 4 for <x a bs = 3 mb for a constant 
octet cross section, independent of rapidity, and color singlet absorption at y — —2 where 
the J/ip can still hadronize in the target at ^/Snn = 200 GeV. The minimum bias results 
are indicated by the dotted lines. Absorption is stronger at central impact parameters but 
dies away gradually at large values of b due to the long tails of the density distributions, 
particularly of deuterium. 

The observed J ftp yield includes feed down from Xcj and tp' decays, giving 

S}%0 - s, z>) = 0.58^ dir (6 - s, z>) + 0.35£(? - s, z') + 0.12^ s (6 - s, z') . (9) 

As discussed previously, in color singlet production, the final state absorption cross section 
depends on the size of the cc pair as it traverses the nucleus, allowing absorption to be 
effective only while the cross section is growing toward its asymptotic size inside the target. 
On the other hand, if the cc is produced as a color octet, hadronization will occur only after 
the pair has traversed the target except at very backward rapidity. We have considered 
a constant octet cross section, as well as one that reverts to a color singlet at backward 
rapidities. For singlets, Sfj^ dir ^ ^ Sffi s since af bs > a*£ s > a' J J^ but, with octets, one 
assumes that Sfj s ^ dir = S^™ = S!^ s since the asymptotic absorption cross section is predicted 
to be equal for all charmonium states. Finally, we have also considered a combination of 
octet and singlet absorption in the context of NRQCD. See Ref. [29] for a detailed description 
of all the absorption models discussed. 

The left-hand side of Fig. 5 shows EKS98 shadowing combined with the absorption 
models described in the text in minimum bias d+Au collisions at ^Snn = 200 GeV. The 
difference between the constant and growing octet cross sections is quite small at large y/Sfm 
with only a minor octet-to-singlet conversion effect at y < —2. Color singlet absorption is 
also important only at similar negative rapidities and is otherwise not different from shad- 
owing alone. The relative combination of nonperturbative singlet and octet contributions in 
NRQCD changes the shape of the shadowing ratio slightly. The large Xcj singlet feeddown 
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contribution reduces the overall absorption effect on inclusive J/ip production. Octet to 
singlet conversion, together with singlet absorption, can be observed at y < — 2 as in the 
other cases. 

Several values of the asymptotic absorption cross section, a a b s = 1, 3 and 5 mb, corre- 
sponding to a = 0.98, 0.95 and 0.92 respectively using Eq. (9) are shown in Fig. 5. These 
values of <7 a b s are somewhat smaller than those obtained for the sharp sphere approximation 
in Eq. (8). The diffuse surface of a real nucleus and the longer range of the density distri- 
bution results in a smaller value of a a b s than that found for a sharp sphere nucleus. There 
is good agreement with the trend of the preliminary PHENIX data [30] for cr abs = — 3 mb. 
We use a value of 3 mb in our further calculations to illustrate the relative importance of 
absorption and shadowing. 

The right-hand side of Fig. 5 compares the EKS98 and FGS parameterizations, all with 
a growing octet cross section with an asymptotic value of <7 a b s = 3 mb. In the region that 
PHENIX can measure, the EKS98 and FGS1 results are essentially indistinguishable. The 
FGSh result lies between the FGSo and EKS98 results at forward rapidity but is also quite 
similar to EKS98 at negative rapidity. 

Figure 6 shows the same calculations for d+Pb collisions at y/S^N — 6.2 TeV at the 
LHC. Although we have plotted the results relative to pp at the same energy, there are 
currently no plans to run LHC pp collisions at energies lower than 14 TeV. Short, lower 
energy proton runs might better establish the energy excitation functions for other processes 
so this possibility is not excluded. In any case, there should be sufficient 14 TeV data to 
produce a high statistics J/ip baseline at this energy. Between the Tevatron Run II data at 
1.96 TeV and the forthcoming LHC pp data at 14 TeV, it should be possible extrapolate to 
6.2 TeV with relative accuracy. 

In addition, it will not be possible to directly compare these results to LHC data unless 
the entire J/ip pt spectrum can be measured since these calculations are integrated over p?. 
Tevatron Run II has shown that it is possible to measure the entire J / ip pt distribution at 
collider energies. A minimum p T cut would reduce the amount of shadowing due to the Q 2 
dependence. 

At 6.2 TeV, the difference between the constant and growing octet cross sections is negli- 
gible for all rapidities shown while only a small absorption effect is seen for the color singlet 
model at y < —5. On the other hand, the difference between the shadowing parameteri- 
zations is larger. Indeed, for y > 2 there is nearly a factor of two difference between the 
EKS98 and FGS ratios, 0.5 and 0.25 respectively. This difference should be measurable in 
the ALICE muon arm, covering the region 2.5 < y < 4. The results at negative rapidity are 
harder to discriminate since the difference between the shadowing ratios is decreased. As 
the antishadowing region is approached, around y = —5, the difference is almost negligible. 

In central collisions, inhomogeneous shadowing is stronger than the homogeneous re- 
sult. The stronger the homogeneous shadowing, the larger the inhomogeneity. In peripheral 
collisions, inhomogeneous effects are weaker than the homogeneous results but some shad- 
owing is still present. Shadowing persists in part because the density in a heavy nucleus is 
large and approximately constant except close to the surface and also because the deuteron 
wave function has a long tail. We also expect absorption to be a stronger effect in cen- 
tral collisions. In Fig. 7, we show the inhomogeneous shadowing and absorption results for 
EKS98 and a abs = 3 mb at v / 3jv7v = 200 GeV as a function of bj Ra for the inhomogeneous 
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to minimum bias ratio (dAu(6)/pp)/(dAu(ave)/pp) =dAu(6)/dAu(ave), shown on the left- 
hand side. The ratios are shown for several values of rapidity to represent the behavior in 
the anti-shadowing (large negative y), shadowing (large positive y) and transition regions 
(midrapidity). The ratios are all less than unity for b/RA < 0.7, with stronger than aver- 
age shadowing and absorption, and rise above unity for large b/RA, weaker than average 
shadowing and absorption. 

We note that for both color octet and color singlet production and absorption, the 
spatial dependence of the absorption, shown in Fig. 4, is weaker than the effect due to 
inhomogeneous shadowing at x ~ 5 x 10~ 4 in Fig. 3. The inhomogeneous shadowing effect 
has a stronger y dependence than absorption which, at collider energies, as illustrated in 
Figs. 5 and 6, is independent of y except at large, negative rapidity. This y-independence 
can help disentangle the effects of shadowing and absorption as a function of centrality. 
For example, if the b dependence of the dAu/pp ratio at RHIC could be determined at 
y ~ —0.5, where the result with <r abs = 0, shown in Fig. 5, indicates dAu/pp ~ 1, the effect 
of absorption alone could be fixed for all measurable y, allowing the effect of shadowing 
alone to be determined for other values of y. 

The right-hand side of Fig. 7 shows the dAu/pp ratios for the same rapidity values as a 
function of the number of binary iViV collisions, iV co n, 



where Ta and Tb are the nuclear thickness functions and a'§ N is the inelastic nucleon-nucleon 
cross section, 42 mb at RHIC. 

The dependence of the ratios on iV co n is almost linear. We do not show results for 
Aeon < 1, corresponding to b/RA > 1.3 on the left-hand side, the point where the 
dAu(6)/dAu(ave) ratios begin to flatten out. The weakest N co n dependence occurs in the 
antishadowing region, illustrated by the y = — 2 result (dot-dashed curve). The trends of the 
ratios as a function of 7V co n are consistent with the PHENIX data from the north muon arm 
(y = 2) and the electron arms (y = 0) but the preliminary PHENIX results from the south 
arm (y = —2) are much stronger than our predictions and, in fact, go the opposite way. The 
overall dependence on N co n is stronger than that obtained from shadowing alone, described 
in Ref. [5] where inhomogeneous shadowing effects depend strongly on the amount of ho- 
mogeneous shadowing. Relatively large effects at low x are accompanied by the strongest b 
dependence. In the transition region around midrapidity at RHIC, the b dependence of the 
ratio dAu/pp due to shadowing is nearly negligible and almost all the N co n dependence at 
y ~ can be attributed to absorption. The effect of absorption alone as a function of iV co ii 
is shown on the right-hand side of Fig. 4. The y = —2 results for color singlet production 
and absorption, in the antishadowing region, are fairly independent of N co \\. 

Figure 8 compares the EKS98 dAu/pp ratios as a function of iV coll to results with the 
FGS parameterizations for the same rapidities as in Fig. 7. Since EKS98 has the weakest 
low x shadowing, given the previous discussion it is not surprising that it also has the 
weakest dependence on iV co ii. The FGSo results have the strongest iV co ii dependence due to 
its strongest overall shadowing. The FGSo ratio goes above unity at y = — 2 since, even 
with a 3 mb absorption cross section, the minimum bias ratio is still greater than one at this 
rapidity. While the FGSh and FGS1 minimum bias ratios are intermediate to the EKS98 
and FGSo ratios, they are more similar to EKS98 at forward rapidities so that, especially for 
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FGS1, the central results are quite similar to those of EKS98 at all rapidities. However, as b 
increases and iV co n decreases, they begin to differ. Since the FGSh and FGS1 inhomogeneous 
parameterizations are tuned to become unity at fO fm, the shadowing component vanishes 
as N co n — > f . On the other hand, there is some residual absorption due to the overlapping 
tails of the deuteron and nuclear density distributions. Thus the y = and y = 2 curves 
meet at N co \\ = 1. At still higher impact parameters, where A^ co u < 1, the y = —2 curve 
intersects the higher rapidity curves. 

While shadowing alone is not incompatible with the preliminary minimum bias PHENIX 
data [30], some absorption should be present. The 3 mb color octet absorption cross section 
is also compatible with the minimum bias data [31] and give reasonable agreement with 
the E866 data when shadowing and other effects, important at forward xf, are included 
[13]. After the data is binned in centrality, some discrimination between models may be 
possible. Including inhomogeneous absorption steepens the dependence on iV co n more than 
with shadowing alone [5]. The preliminary PHENIX centrality dependence [30] seems to 
be relatively flat for the central electron arms and the north muon arm (positive rapidity) 
while the centrality dependence in the south muon arm (negative rapidity) seems to increase 
strongly with iV co ii. The EKS98 and FGS1 calculations are most compatible with the central 
and north muon arm data while FGSo and FGSh seem somewhat too strong relative to the 
data. Perhaps with a smaller absorption cross section, these calculations could be made 
more compatible with the data. None of the calculations agree with the preliminary south 
muon arm data, exhibiting very little dependence on iV co ii. The final PHENIX data as 
well as any data from a future d+Au run would help better determine both the amount of 
gluon shadowing as well as the relative contributions of shadowing and absorption, both in 
minimum bias collisions and as a function of centrality. It may be that the PHENIX data 
would be more compatible with color singlet rather than color octet absorption although 
this would contradict fixed target results where shadowing is weak [1,10]. 

In Figs. 9 and 10 we present our inhomogeneous shadowing and absorption calculations 
for d+Pb collisions at \/Snn — 6.2 TeV at the LHC. Figure 9 shows only the EKS98 results 
while Fig. 10 compares the EKS98 results as a function of N co n to those of FGS. Here we 
have also included results for y = ±4, within the range of the ALICE muon arm. Given 
that the rapidity range of the muon arm encompasses the crossover point where dPb/pp ~ 1 
at y ~ —3.9, the centrality dependence of absorption alone could be determined and used 
to calibrate the inhomogeneous shadowing effects. Note that is is only possible to reach 
y ~ —3.9 in ALICE by switching the beam directions and running Pb+d since the muon 
arm is only on one side of midrapidity. Both ALICE and CMS should be able to measure 
J/ijj production at y — ±2 and 0. Since the exact pp baseline is likely to be absent, as 
discussed above, perhaps a better baseline for the centrality dependence would be the ratio 
relative to the minimum bias result, shown on the left-hand side of Fig. 9 as a function 
of b/RA- The change in dPb(6)/dPb(ave) with b/RA at the LHC is stronger than that of 
dAu(6) /dAu(ave) at RHIC, shown in Fig. 7. 

The results for y — 4, ±2 and 0, all in the low x shadowing region, are rather closely 
grouped together. This should not be surprising because the EKS98 shadowing ratios shown 
in Fig. 6 are not very strong functions of rapidity. The grouping of the results is even more 
striking as a function of iV co n, shown on the right-hand side of Fig. 9. Note also here that the 
x-axis scale is expanded relative to that of RHIC since a^ N = 76 mb at 6.2 TeV, increasing 
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the number of inelastic NN collisions possible at the LHC The expanded A^ co ii scale also 
has implications for the FGSh and FGS1 parameterizations with their sharp 10 fm cutoff 
on shadowing. As seen in Fig. 10(c) and (d), the curvature of these parameterizations is 
much stronger than the linear dependence produced by Eq. (6). This increased strength 
is due to the larger value of ct^n 

at the LHC. Now iV coU « 1 at b ~ 10 fm. Here FGSh 
and FGS1 shadowing vanishes and the curves all come together at the value of the ratio 
corresponding to the amount of residual absorption. The curvature increases with y since 
the larger shadowing ratio has to increase to unity faster than those at negative rapidity 
where shadowing is weaker. 

The results in Figs. 9 and 10 show that while, at y < in Fig. 6 good statistical accuracy is 
necessary to discriminate between shadowing models, the difference in their 7V co n dependence 
may be sufficient to do so. 

Finally, since it is not yet clear whether proton or deuterium comparison runs will be 
used at the LHC, Fig. 11 shows the results for p+Pb collisions at y/S^N = 8.8 TeV, the 
appropriate energy for proton-nucleus collisions. We also have plotted the results assuming 
that the center of mass rapidity is at y — 0, as is the case for symmetric beams. However, 
since the charge-to-mass, Z/A, ratio is different for protons and lead, the proton beam 
has an energy of 7 TeV while the Pb beam has energy 2.75 TeV. The asymmetric energies 
lead to a shift of the center of mass rapidity of 0.5 units in the direction of the proton 
beam. Since the difference between the Z/A ratios of deuterium and lead is much smaller, 
the shift in d+Pb collisions is only 0.1 unit. The p+Pb antishadowing peak is shifted 
further toward negative rapidity while the region where the parameterizations flatten begins 
closer to midrapidity than in Fig. 6. Thus the EKS98 results as function of N co \\, shown 
on the right-hand side of Fig. 11 for y = 4 (solid curve) and y = (dashed curve) are 
quite close together: the shadowing ratio changes very little as a function of rapidity. The 
corresponding FGSh results are shown in the dot-dashed and dotted curves respectively. In 
both cases, the iV co ii dependence is stronger than in d+Pb collisions, as exhibited by the 
disappearance of both shadowing and absorption as N co n — > 1 for the two inhomogeneous 
shadowing parameterizations. The change in slope is particularly strong for EKS98 because 
the proton is treated as a point particle while the deuteron is an extended system. The 
maximum iV co ii for p+Pb is slightly larger than half that for d+Pb. If not for the higher 
inelastic cross section, u^n = 81 mb at 8.8 TeV, iV coll in p+Pb would be exactly half that 
of d+Pb. 

In conclusion, the preliminary PHENIX data show that J / ip production in d+Au collision 
is modified with respect to production in pp collisions at the same energy. This modification 
is consistent with initial-state shadowing plus final-state absorption seen in fixed-target 
experiments at 800 GeV [13]. More precise measurements may help to better set the level of 
absorption allowed by the data. Precision measurements of the centrality dependence may 
also help. Corresponding data from the LHC will be more useful for separating and refining 
shadowing models due the very low x range available for large rapidity measurements. 

We thank K.J. Eskola and V. Guzey for their shadowing parameterizations. We also 
thank S.R. Klein and M.J. Leitch for discussions. This work was supported in part by the 
Division of Nuclear Physics of the Office of High Energy and Nuclear Physics of the U.S. 
Department of Energy under Contract Number DE-AC03-76SF00098. 
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FIGURES 
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FIG. 1. The shadowing parameterizations are compared at the scale \x = 2m c = 2.4 GeV. The 
solid curves are EKS98, the dashed, FGSo, dot-dashed, FGSh, and dotted, FGS1. 
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FIG. 2. We give the average value of the nucleoli momentum fraction, X2, of J/ip production 
in pp collisions as a function of rapidity for (a) the CERN SPS with y/S = 19.4 GeV, (b) RHIC 
with y/S = 200 GeV and (c) the LHC with s/S = 6.2 TeV. 
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FIG. 3. The WS (dot-dashed) and p (solid) inhomogeneous shadowing parameterizations are 
compared to the inhomogeneous FGS shadowing parameterization (dashed) at the same value of 
the homogeneous ratio, indicated by the horizontal solid line. 
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FIG. 4. The J/^> dAu/pp ratio for absorption alone with <7 a b s = 3 mb as a function of impact 
parameter (left-hand side) and as a function of the number of nucleon-nucleon collisions, N co \\, 
(right-hand side) for a constant octet (all y), solid, and singlet (y = — 2), dashed. The homogeneous 
results are indicated by the dotted lines. 
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FIG. 5. Left-hand side: The J/ip dAu/pp ratio with EKS98 at 200 GeV as a function of rapidity 

for (a) constant octet, (b) growing octet, (c) singlet, all calculated in the CEM and (d) NRQCD. For 

(a)-(c), the curves are no absorption (solid), a a ^ s = 1 (dashed), 3 (dot-dashed) and 5 mb (dotted). 

For (d), we show no absorption (solid), 1 mb octet/ 1 mb singlet (dashed), 3 mb octet/3 mb singlet 

(dot-dashed), and 5 mb octet/3 mb singlet (dotted). Right-hand side: The J/ip dAu/pp ratio at 

200 GeV for a growing octet with a a b s = 3 mb is compared for four shadowing parameterizations. 

We show the EKS98 (solid), FGSo (dashed), FGSh (dot-dashed) and FGS1 (dotted) results as a 

function of rapidity. 
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FIG. 6. Left-hand side: The J/ip dPb/pp ratio with EKS98 at 6.2 TeV as a function of rapidity 

for (a) constant octet, (b) growing octet, (c) singlet, all calculated in the CEM and (d) NRQCD. For 

(a)-(c), the curves are no absorption (solid), a a ^ s = 1 (dashed), 3 (dot-dashed) and 5 mb (dotted). 

For (d), we show no absorption (solid), 1 mb octet/ 1 mb singlet (dashed), 3 mb octet/3 mb singlet 

(dot-dashed), and 5 mb octet/3 mb singlet (dotted). Right-hand side: The J/ip dAu/pp ratio at 

200 GeV for a growing octet with a a b s = 3 mb is compared for four shadowing parameterizations. 

We show the EKS98 (solid), FGSo (dashed), FGSh (dot-dashed) and FGS1 (dotted) results as a 

function of rapidity. 
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FIG. 7. Left-hand side: The J/^> ratio dAu(6)/dAu(ave) as a function of b/RA- Right-hand 
side: The ratio dAu/pp as a function of N co \\. Results are shown for y = —2 (dot-dashed), y = 
(dashed) and y = 2 (solid) at 200 GeV for a growing octet with o" a b s = 3 mb and the EKS98 
parameterization. 
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FIG. 8. The ratio dAu/pp as a function of 7V coU for the EKS98 (a), FGSo (b), FGSh (c) 
and FGS1 (d) shadowing parameterizations. The calculations with EKS98 and FGSo use the 
inhomogeneous path length parameterization while that obtained by FGS is used with FGSh and 
FGS1. Results are given for y = —2 (dot-dashed), y = (dashed) and y = 2 (solid) at 200 GeV for 
a growing octet with <7 a b s = 3 mb. 
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FIG. 9. Left-hand side: The J/tp ratio dPb(6)/dPb(ave) as a function of b/RA- Right-hand 
side: The ratio dPb/pp as a function of N co \\. Results are shown for y = —4 (dot-dot-dot-dashed), 
y = —2 (dotted), y = (dot-dashed), y = 2 (dashed) and y = 4 (solid) at 6.2 TeV for a growing 
octet with o" a bs = 3 mb and the EKS98 parameterization. 
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FIG. 10. The ratio dPb jpp as a function of 7V coll for the EKS98 (a), FGSo (b), FGSh (c) 
and FGS1 (d) shadowing parameterizations. The calculations with EKS98 and FGSo use the 
inhomogeneous path length parameterization while that obtained by FGS is used with FGSh and 
FGS1. Results are given for y = —4 (dot-dot-dot-dashed), y = — 2 (dotted), y = (dot-dashed), 
y = 2 (dashed) and y = 4 (solid) at 6.2 TeV for a growing octet with cr a bs = 3 mb. 
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FIG. 11. Left-hand side: The ratio pPb/pp at 8.8 TeV as a function of rapidity for the EKS98 
(solid), FGSo (dashed), FGSh (dot-dashed) and FGS1 (dotted) shadowing parameterizations with 
a 3 mb octet absorption cross section. Right-hand side: The N co \\ dependence for EKS98 (y = 4, 
solid curve and y = 0, dashed curve) and FGSh (y = 4, dot-dashed curve and y = 0, dotted curve). 
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